Self-energy of image states on copper surfaces 
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We report extensive calculations of the imaginary part of the electron self-energy in the vicinity of 
the (100) and (111) surfaces of Cu. The quasiparticle self-energy is computed by going beyond a 
free-electron description of the metal surface, either within the GW approximation of many-body 
theory or with inclusion, within the GWF approximation, of short-range exchange-correlation effects. 
Calculations of the decay rate of the first three image states on Cu(lOO) and the first image state 
on Cu(lll) are also reported, and the impact of both band structure and many-body effects on the 
electron relaxation process is discussed. 



I. INTRODUCTION 

It is "well-known that an electron outside a metal sur- 
face experiences an effective potential with the asymp- 
totic image form Vim = — l/4(z — zq), zq representing the 
image-plane position. If the bulk band structure pro- 
jected into the surface presents a band gap near the vac- 
uum level {Ey = 0), an electron with energy E < can be 
trapped in the potential-well created on the vacuum side 
of the surface by the gap and the image potential. The 
long range character of the image potential gives rise to a 
series of unoccupied Rydbftrg-like states, which converge 
towards the vacuum leveLEJa These[Sa-called image states 
were first detected experimentallyuu by the techniques 
of inverse photoemission,Ellj and the first high-resolution 
measurements of image states were, made by the use of 
two-photon photoemission (2PPE).Lriij 

Image-potential states are localized in the vacuum re- 
gion of the surface. Hence, they are almost decoupled 
from bulk states and are expected to have much longer 
lifetimes than bulk excitations. Lifetimes of image states 
had been determined from their spectral width in 2PPE 
spectroscopy, and the advent of time-resolved 2PPE (TR- 
TPPE) has enabled a direct measurement of image-jSta^ 
lifetimes on-the, (100) and (111) surfaces of silverll3~tll 
and copper JIjIIj 

Calculations of image-state lifetimes were reported in 
Refs. ^ and |2l|, within a many-body free-electron de- 
scription of the metal surface and with the use of simpli- 
fied models to approximate both initial and final elec- 
tronic states and, also, the screened Coulomb inter- 
action. Later on, the role that crystal-induced sur- 
face states, not present within a free-electron descrip- 
tion of the spLLd. play in the decay of image states was 
invcstigated.t3c3 Self-consistent many-body calculations 
of image-state lifetimes that go beyond a free-electron de- 
scription of the metal surface have beea. carried out only 
very recently on coppercil and lithiumcS surfaces. 

In this paper we report extensive calculations of the 



imaginary part of the image-electron self-energy in the 
vicinity of the (100) and (111) surfaces of copper, which 
we conipute within the GW approximation of many-body 
theory.tj Single-particle wave functions are obtained by 
solving the Schrodinger equation with a realistic one- 
dimensional model potential, and the screened interac- 
tion is evaluated within the random-phase approxima- 
tion (RPA).EZl Then, we present the results of calcula- 
tions of the lifetime of the first three image states on 
Cu(lOO) and the first image state on Cu(lll). These 
calculations are carried out within a GW-RPA descrip- 
tion of the self-energy and, —also, with use of the so- 
called GWr approximation,E3'E2l which includes short- 
range exchange-correlation (XC) effects not present in 
the GW-RPA. We also consider simplified models for 
the description of both the electronic wave functions and 
the screened interaction, we discuss the impact of band 
structure and many-body effects on image-state lifetimes, 
and account for the potential variation parallel to the 
surface through the introduction of the effective mass. 
We present systematic investigations of the role that all 
quantities involved in the evaluation of the electron self- 
energy play in the coupling of image states with the solid. 
We find that a detailed description of thes quantities is 
of crucial importance in the understanding of the origin 
and magnitude of decay rates of image states. Finally, 
differences between decay rates of image states on the 
(100) and (111) surfaces of Cu are investigated, and we 
focus our attention on the role that the various available 
bulk and surface states play in the electron relaxation 
process. 

Unless otherwise is stated, atomic units are used 
throughout, i. e., ~ fi ^ rric = 1. 

II. THEORY 

Take an inhomogeneous electron system. The damping 
rate of an excited electron in the state V'o(r) with energy 
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-Bo is obtained as the projection of the imaginary part of 
the electron self-energy, E(r,r',£^o), over the state itself: 

= -2 J dr J dr'V'5(r)ImE(r,r';£;o) V^o(r')- (1) 

We consider a periodic-film model of the solid. The 
film is taken to be translationally invariant in the plane 
of the surface, which is assumed to be normal to the 
z axis, and departure of motion along the surface from 
free-electron behaviour is accounted through the intro- 
duction of the effective mass. Hence, single-particle wave 
functions are taken to be of the form 



with energies 



^11 

2 mo 



(2) 



(3) 



The wave function (/>o(^) and energy Eq describe motion 
normal to the surface, ky is a wave vector parallel to 
the surface, A is the normalization area, and mo is the 
effective mass. Introduction of Eq. (||) into Eq. (|l|) yields 
the following expression for the damping rate: 

r"^ = -2 J dz J dz (t)Q{z)linJ:{z,z;kii,Eo)(j)oiz'), 

(4) 

where E(z, z'; k|| , i?o) represents the two-dimensional 
Fourier transform of the electron self-energy. 

In the so-called GW approximation, one considers only 
the first-order term in a series expansion of the self- 
energy in terms of the screened interaction, 

^^2^7 (2^ ^(^'^'^11'^) 
X iG(z,z';k|| -q||,£:o -^^), (5) 

and after replacing the Green function G(z, z'; k|| , a;) by 
that of non-interacting electrons, one finds 



ImS(z,z';k||,i?o) = 5]0}(z')0/(z) 
/ 
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The sum in Eq. is extended over a complete set 
of final states 0/(z) with energy e/ {Ep < Ef < Eq), 
Ep is the Fermi energy, and VF(z, z'; qy , is the two- 
dimensional Fourier transform of the screened interac- 
tion: 



VK(z,z';q||,w) = w(z,z';q||) + J dzi J dz2 

xv{z, zi; q||) x(2:i, 22; q||, t^) v{z2,z'; q\\), (8) 

with w(z,z'; q||) and xiz, z';q_\^,uj) representing the bare 
Coulomb interaction, 



i;(z,z';q||) = ^e-««l^-^'l. 



(9) 



and the density-response function, respectively. 

Within RPA, x(z, z'; qy , oj) satisfies the integral equa- 
tion 



X{z,z';q\\,uj) = xo{z,z',q\\,uj) + J dzi J dzi 

xxo(^, zi, q|| , w) u(zi, Z2; q|| ) x(^2, z', q|| , t^), (10) 

where xq{^z,z' ,<\\\,ld) represents the density-response 
function of non-interacting electrons,tj 

Xo(2:,2:',q||,w) = ^0((z)0;',(z)0;/(z')(/'r(z') 
l,V 

dk|| Q(Ef-Ei)-(){Ep-Ei,) 



with 



and 
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and Q{x) being the Heaviside step function. 

In the so-called GWF approximation, which includes 
XC effects not present in the GW-RPA, the electron self- 
energy is of the GW form, i. c., it is given by Eq. (H), 
but with an effective screened interaction 



W^(z,z';q||,tj) = w(z,z';q||) + j dzi y dza 
X [w(zi,Z2;q||) K.j,c{zx,Z2;<\\\,lS)\ 
xx(2i, Z2; q||, w) w(z2, z'; q||), (14) 

where 

X(z,z';q||,w) = xo(^,^',q||,'^) + ^ dzi ^ dz2 
XXo(2:,2:i,q||,w) [w(z,zi;q||) 

+K^J^zi,Z2\(\\\,'jj)\ x(22,z',q||,tj). (15) 

Here, K^^J^z^z' represents the two-dimensional 
Fourier transform of the XC kernel K^(iv,v' \u}), which 
accounts through Eqs. (^ and ( |l^ ) for the reduction in 
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the electron-electron interaction due to the existence of 
short-range XC effects associated to the probe electron 
and to screening electrons, respectively. In tlie static 
liniit-(tj 0), density- functional theory (DFT)lII shows 

thatka 



dn{r)dn{r') 



(16) 



"o(r) 



where i5a;c[n] represents the XC energy functional and 
no(r) is the actual density of the electron system. In 
the local-density approximation (LDA), the XC kernel is 
approximated by a contact S function, and one finds 



K^ciz,z';q\i,uj ^ 0) = 



dn^ 



S{z~z'), (17) 



no{z) 



where E^dn] now represents the SC energy of a homo- 
geneous electron gas of density nB Introduction of this 
static XC kernel into Eq. (|l^) represents an adiabatic 
extension of the LDA to finite frequencies j^and yields the 
so-called time-dependent LDA (TDLDA).bl 

The single-particle wave functions (t>i{z) entering Eqs. 
®; ® (01) ^^"^ simply eigenfunctions of a one- 
electron hamiltonian. In the particular case of the 
RPA, 4>i{z) are self-consistent eigenfunctions of the one- 
electron Hartree hamiltonian, and within TDLDA they 
are obtained by solving the Kohn-Sham equation of DFT 
with use of the LDA XC potential. 



V^z) 



dn 



(18) 



no(z) 



Neither the Hartree self-consistent eigenfunctions nor 
the LDA wave functions produce the correct image-like 
asymptotic potential behaviour on the vacuum side of 
the surface. 

For a realistic description of the metal surface, we solve 
for (j)i{z) a single-particle time-independent Schrodinger 
equation with the one-dimensional model potential sug- 
gested in Ref. which approaches, far outside the sur- 
face, the classical image potential. This one-dimensional 
potential has four adjustable parameters, which are cho- 
sen so as to reproduce the width and position of the en- 
ergy gap at the T point (k|| = 0) and, also, the bind- 
ing energies of both the n — crystal-induced surface 
state at F and the first {n — I) image-potential induced 
state. Probability amplitudes of the n — 1 image state on 
Li(llO) and Cu(lOO), as obtained with use of this model 
potential, have been found to be in|-e}ccellent agreemept 
with first principles pseudopotentialEJ and all-electronE3 
calculations, respectively. 



III. RESULTS AND DISCUSSION 

Input of our calculation of image-state decay rates, as 
jiven by Eq. (0), are the image-state wave function (j)o{z) 



and the image-electron self-energy S(z, z'; k|| , _Eo). We 
use films of up to 50 layers of atoms and 80 interlayer- 
spacing vacuum intervals, thereby ensuring that finite- 
slab effects are negligible. Image-state wave functions are 
taken to be eigenfunctions of the one-dimensional model 
hamiltonian described above (MP). For the evaluation of 
the image-electron self-energy, we use in Eq. (|^) either 
the MP wave functions or the self-consistent jellium LDA 
eigenfunctions of the one-electron Kohn-Sham hamilto- 
nian without (J) and with (JG) the restriction that only 
final states with energy e / lying below the projected band 
gap are allowed. The screened interaction entering Eq. 
(P) is evaluated either within tha-specular-reflexion model 
(SRM) of Ritchie and Marusak,EZl with use of the approx- 
imate-surface response function of Persson and Zaremba 
(PZ)J23 or within the self-consistent approaches of Eqs. 
(^ and (|l^ ) with the single-particle eigenstates entering 
Eq. ( prj ) being MP wave functions. 

Probability densities, |0o(-z)P, for the n — 1 image 
state on the (100) and (111) surfaces of Cu, as obtained 
with use of the one-dimensional model potential of Ref. 
|35| , are represented in Fig. 1. In the case of Cu(lOO) 
(dotted line), the n = 1 probability-density has a maxi- 
mum at 3.8 A outside the crystal edge (z = 0), which we 
choose to be located half a lattice spacing beyond the last 
atomic layer, and the penetration into the bulk (z < 0) 
crystal is found to be of 5%. On the (100) surface of 
Cu, the n = 1 imasc, state is close to the center of the 
projected band gap,E3 which results in a very small pen- 
etration into the bulk. However, the n = 1 image state 
on Gii(lll) (solid line) is located right at the top of the 
gapjlj the solution in the bulk is an s-like wave function 
with the matching at the surface occurring at minimum 
amplitude, the maximum of the probability-density in 
the vacuum is closer from the surface, at 2.3 A, than in 
the case of Cu(lOO), and the penetration into the bulk 
is found to be of 22%. The first image state on Cu(lOO) 
and Cu(lll) has binding energies of 0.57 and 0.82 eV, 
respectively. 

Fig. 2 shows full GW-RPA calculations of the imag- 
inary part of the ?i = 1 image-electron self-energy 
Im [— ^(z, z'; kj[ — 0, Eq)] in the vicinity of the (100) and 
(111) surfaces of Cu, as obtained from Eq. ^ with use 
of MP wave functions both in Eq. (^ and Eq. (|ll|) and 
with all effective masses set equal to the free-electron 
mass. The imaginary part of the electron self-energy is 
represented in this figure as a function of z and for a fixed 
value of z' . In the top panel, z' is fixed at about three 
atomic layers {z' ^ —5 A) within the bulk, showing that 
Im(— E) has a maximum at z = z', as in the case of a 
homogeneous electron gas. When z' is fixed at the crys- 
tal edge (z' '--^ 0), as shown in the middle panel of Fig. 2, 
we find that Im(— S) is still maximum at z = z', but the 
magnitude of this maximum now being enhanced. The 
bottom panel of Fig. 2 corresponds to z' being fixed at 
about three atomic layers (z' ~ 5A) from the surface 
into the vacuum. In this case the maximum magnitude 
of Im(— E) occurs at z ~ rather than for z — z', show- 
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ing a highly nonlocal behaviour of the imaginary part of 
the electron self-energy at the surface. This nonlocalifoj 
of the cusp of Im(— E) was also shown by Deisz et a/,E3 
within a free-electron (jellium) description of the surface. 

The magnitude of the maximum of Im(— E) is plotted 
in Fig. 3, as a function of z', showing that it is an oscillat- 
ing function of z' within the bulk and reaches its highest 
value near the surface. The oscillatory behaviour within 
the bulk is dictated by the periodicity of final-state wave 
functions 4>f{z) in the crystal, and the highest value near 
the crystal edge is the result of electron-hole pair cre- 
ation taking place mainly in the vicinity of the surface. 
We note that the magnitude of Im(— S) near the surface 
is larger for Cu(lll) than for Cu(lOO), although the band 
gap on Cu(lll) extends below the Fermi level, thus the 
available phase space on this surface becoming highly re- 
stricted. However, while the crystal-induced n = sur- 
face state on Cu(lOO) is located outside the projected 
band gap and represents, therefore, a so-called surface 
resonance, the n — Q surface state on Cu(lll) is located 
within the band gap and provides an important channel 
for the decay of image states. As a result, the imaginary 
part of the image-electron self-energy near the (111) sur- 
face of Cu is largely enhanced. This is illustrated in Fig. 
4, where contributions to the magnitude of the maximum 
Im(— E) of the n = 1 image state on Cu(lll) have been 
plotted separately, according to whether transitions to 
bulk states (dotted line) or to the n — surface state 
(dashed line) occur. 

Now we focus on the evaluation of image-state life- 
times, and we set the wave vector of the image electron 
parallel to the surface, k|| , equal to zero. Coupling of the 
image state with the crystal occurs through the penetra- 
tion of the image-state wave function into the solid and, 
also, through the evanescent tails of bulk states outside 
the crystal. Accordingly, we have calculated separately 
the various contributions to the damping rate by con- 
fining the integral in Eq. (^ to either bulk {z < 0) or 
vacuum (z > 0) coordinates: 



'bulk 



(19) 



where rj~}j,, t^Jc ^'^d t^]^^ represent bulk, vacuum and 
interference contributions, respectively. The results of 
our calculations for the decay rate of the n — 1 image 
state on the (100) and (111) surfaces of Cu are presented 
in Tables I, II, and III. 

In order to investigate the impact of band structure 
effects on the damping rate of image states, we present 
in Table I the results of our full GW-RPA calculations 
for the decay rate of the n — 1 image state on Cu(lOO), 
as obtained from Eq. (0) with use of either J, JG or MP 
wave functions in Eq. (ph, with use of MP wave fuctions 
in Eq. (pi]), and with all effective masses set equal to 
the free-electron mass. For g|| > y/2{Eo — Eg) {Eg rep- 
resents the bottom of the projected band gap) all final 
states with energy Ef < Eq lie below the gap, thereby 
bulk and interference contributions to the decay rate re- 



maining nearly unaffected by this restriction. However, 
as the coupling of the image state with the crystal oc- 
curring through the tails of bulk states outside the crys- 
tal is expected to be dominated by vertical transitions 
(Q|| ^ 0)i the vacuum contribution to the decay rate be- 
comes noticeably smaller as final states lying within the 
projected band gap are not allowed, and this restriction 
results in a total decay rate that is smaller by a factor of 
~ 3. However, the decay rate is nearly insensitive to the 
actual choice of the one-particle wave functions entering 
Eq. (||), as long as only final states lying below the pro- 
jected band gap are included in the jellium calculation. 
We have also performed calculations with use of either J 
or JG wave functions in Eq. (11), and have found that 
the decay rate is insensitive to the actual choice of the 
one-particle wave functions entering Eq. (0). In the 
case of the (111) surface of Cu, the bottom of the pro- 
jected band gap is located just below the Fermi level, and 
both the existence of the n — surface state, not present 
within a jellium model of the surface, and the restricted 
available phase space play a key role in the determination 
of the n = 1 image-state decay rate. 

In Tables II and III we present the results of our cal- 
culations for the decay rate of the n = 1 image state 
on CuflOO) and (111), respectively, as obtained from 
Eq. (P with use of MP wave functions in both Eqs. 
(^) and (^ij). First of all, we set all effective masses 
equal to the free-electron mass, and focus on the role 
that an accurate description of the screened interaction 
plays in the coupling of image states with the solid. 
Hence, we use three different models for the evaluation 
of W^(z, z'; q|| , w). First, the specular-reflexion model 
(SRM) of Ritchie and MarusakEJ has been considered, 
thereby assuming that bulk electrons are specularly re- 
flected at the surface with no interference between ingo- 
ing and outgoing waves. Secondly, for the vacuum con- 
tribution to the decay rate (z > 0, z' > 0) the surface 
response function suggested by Persson and ZarembacS 
(PZ) has been used. Finally, the screened interaction 
has been evaluated from Eq. (||) by solving the RPA in- 
tegral equation (Eq. (^0|)) for the density-response func- 
tion (GW-RPA). 

We note that simplified jellium models (SRM and PZ) 
for the evaluation of the screened interaction yield un- 
realistic results for the image-state lifetime. Bulk con- 
tributions to the linewidth are approximately well de- 
scribed within the specular reflexion model, small differ- 
ences resulting from an approximate description, within 
this model, of the so-called bregenzung effects. How- 
ever, as within this model quantum-mechanical details of 
the surface are ignored, it fails to describe both vacuum 
and interference contributions to the decay rate. These 
quantum-mechanical details of the surface are approxi- 
mately taken into account within the PZ jellium model, 
thereby resulting in a better approximation for the vac- 
uum contribution to the decay rate than the SRM, but 
within the PZ model one neglects the coupling of the 
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image state with the crystal that occurs through the 
penetration of the image-state wave function into the 
sohd. Discrepancies between vacuum contributions ob- 
tained within this model and our more realistic full RPA 
calculations appear as a result of the jellium model of 
Ref. ^ being accurate only for q\\/qF and uj/Ep << 1 
(qp is the Fermi momentum, i. e., Ep = qp/2). 

Now we look at the impact of short-range XC effects, 
which are well-known to reduce electron-electron interac- 
tions both between the image electron and the electron 
gas and between screening electrons themselves. Hence, 
we still set all effective masses equal to the free-electron 
mass, we introduce MP wave functions into Eqs. and 
(|rT|), and compare the GW-RPA calculations described 
above with the results that we obtain from either Eq. (||) 
(GW) or Eq. dl^) (GWr) by using the density-response 
function of Eq. (|15|) with the LDA XC kernel of Eq. ^) 
(TDLDA). An inspection of the results presented in Ta- 
bles I and II indicates that the existence of XC effects 
between screening electrons enhances the decay proba- 
bility of image states. Nevertheless, this enhancement is 
more than compensated by the large reduction of the de- 
cay rate produced by the presence of a XC hole around 
the image-state electron. Consequently, GW-RPA cal- 
culations produce decay rates that are larger than their 
more realistic GWP-TDLDA counterparts by no more 
than ~ 5% in both (100) and (111) surfaces of Cu. 

Finally, we account for potential variation parallel to 
the surface through the introduction of a realistic ef- 
fective mass. The dispersion i?o(k||) of image states 
has been determined experimentally with the use of iiu 
verse photoemission techniques at off-normal emission,t2l 
showing that the effective mass of image states in both 
Cu(lOO) and Cu(lll) are ~ 1, i. e., the free electron 
mass. Measurements of the dispersion of the n = 
surface resonance/state on the (100) and (111) surfaces 
of Cu havCr^iekled effective masses of 0.50 and 0.42, 
respectively,E2rE3 and for bulk states entering Eq. (^ we 
have chosen to increase the effective mass from our com- 
puted value of m/ = 0.40 and 0.22 at the bottom of the 
gap in Cu(lOO) and Cu(lll), respectively, to mj = 1 at 
the bottom of the valence band. As the effective mass of 
all available final states is smaller than the free-electron 
mass, the n — 1 image-state decay rate of both Cu(lOO) 
and Cu(lll) is found to be about 20 — 25% smaller than 
in the case of free-electron behaviour along the surface 
{nif = 1). On the one hand, there is the effect of the de- 
crease of the available phase space, which is easily found 
to scale as y/fnj. On the other hand, as the effective 
mass decreases the decay from the image state occurs, for 
a given energy transfer, through smaller parallel momen- 
tum transfer, which may result in either enlarged or di- 
minished screened interactions, depending on the magni- 
tude of momentum and energy transfers. Our results also 
indicate that GW-RPA and GWF-TDLDA calculations 
nearly coincide, as in the case of free-electron behaviour 
along the surface, which is a consequence of the compe- 



tition between XC effects associated to the image-state 
electron and to screening electrons themselves. GWF- 
TDLDA calculations, as obtained with use of realistic 
effective masses for the description of final-state wave 
functions, yield decay rates of the n = 1 image state on 
Cu(lOO) and Cu(lll) of 17 and 28.5 meV, respectively, in 
excellenL^greement_with the experimeajfcally determined 
lifetimes^ of 40 ± GBM and 22 ± 3 fsO 

With the aim of investigating the role that the vari- 
ous available bulk and surface states play in the decay 
of image-potential states, now we focus on our full GW- 
RPA calculation of the n = 1 image-state decay rate. Fig. 
5 exhibits tJ^ separate contributions to from all fi- 
nal states lying below the projected band gap in Cu(lOO) 
(curves with circles) and Cu(lll) (curves with squares), 
as obtained with the final-state effective mass set equal to 
the free-electron mass (solid lines) and with use of realis- 
tic values for irif (dashed lines). In the case of Cu(lll), 
there is still a large contribution to r"^ from the decay of 
the n = 1 image state into the crystal-induced n = sur- 
face state, lying within the projected band gap, which ap- 
proximately represents 40% of the total decay rate. Figs. 
5(a), (b), (c) and (d) exhibit bulk, vacuum, interference, 
and total contributions to t^^. As the n = 1 image-state 
wave-function overlap with the bulk is larger in Cu(lll) 
than in the case of Cu(lOO), bulk contributions to tJ^ de- 
cay rates are much larger for Cu(lll) than for Cu(lOO), 
as illustrated in Fig. 5(a). However, the large bulk-state 
overlap in Cu(lll) is partially counterbalanced by the 
band gap extending on the (111) surface of Cu below the 
Fermi level. On the other hand, vacuum and interfer- 
ence contributions to the decay rate in both surfaces of 
Cu are comparable in magnitude and opposite in sign, 
yielding total decay rates that differ little from the bulk 
contribution. We also note that the total decay rate in 
Cu(lOO) (22 and 17.5 meV, with to/ = 1 and to/ ^ 1, 
respectively) nearly coincides with the X)/ ''"/^ contribu- 
tion from all bulk states in Cu(lll) (21 and 17nieV). 
Hence, differences between total decay rates in Cu(lOO) 
and Cu(lll) appear as a consequence of the presence of 
the 71 = surface state on Cu(lll), which provides a key 
channel for the decay mechanism. 

Decay rates of image-potential states on Cu(lOO) with 
quantum number n < 3 are presented in Table IV, to- 
gether with the experimentally determined lifetimes re- 
ported in Refs. |l^ and |l^. As before, the wave vector of 
the image state parallel to the surface, k, has been set 
equal to zero, all the wave functions entering Eqs. (|^), 
(^) and ( prl ) have been chosen to be MP wave functions, 
and the screened interaction has been evaluated within 
the GW-RPA. Decay rates of the n = 2 and n = 3 image 
states have been split as follows: 

^"'=E^7'+^n=l (20) 
/ 

and 
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r-'=Y.^j'+r-l,+T-l,, (21) 
/ 

respectively, where tJ^ represents, as in the case of 
the n = 1 image state, the decay rate from transitions to 
all available states lying below the projected band gap, 
and represents the contribution from decay into the 
lower lying n image state. We observe that lower lying 
image states can give noticeable contributions to the de- 
cay rate of excited, i. e., n = 2 and n = 3 image states. 
Decay of these states into the n = 1 image state results in 
linewidths that represent ~ 10% of the total linewidth. 
Decay of the n = 3 image state into the n = 2 lower 
lying state results in a linewidth that represents ^ 3% of 
the total linewidth. We also observe that both our calcu- 
lated and the experimentally determined lifetimes of im- 
age states in Cu(lOO) increase with rt as ~ in agree- 
ment with previous theoretical predictions .u Discrepan- 
cies between our calculated inelastic lifetimes of excited 
image states {n > 1) and experimental measurements 
may be attributed to scattering with phonons, which oc- 
curs on a time scale (~ Ips) that is comparable for these 
states to the electron-electron relaxation time. 



IV. SUMMARY 

We have reported extensive self-consistent calculations 
of the imaginary part of the electron self-energy and the 
decay rate of image states on the (100) and (111) sur- 
faces of Cu, by going beyond a free-electron description 
of the metal surface. We have found that the imagi- 
nary part of the electron self-energy outside the surface 
is highly nonlocal, and have found the magnitude of the 
maximum of this quantity to reach its highest value near 
the surface. We have presented the results of calculations 
of the lifetime of the first three image states on Cu(lOO) 
and the first image state on Cu(lll), and have focused 
on the impact of band-structure and many-body effects 
on these quantities. We have found that band-structure 
effects on the evaluation of final-state wave functions 
may be approximately accounted through introduction 
of the restriction that only final states lying below the 
projected band gap are allowed, and the impact of the 
band structure through the evaluation of the screened 
interaction has been found to be large. We have shown 
that simplified jellium models for the electronic response 
yield unrealistic results for the lifetime of image states, 
thereby a detailed description of the screened interaction 
being of crucial importance in the understanding of the 
origin of image-state lifetimes. We have evaluated the 
screened interaction within three different self-consistent 
many-body schemes, depending on whether XC effects 
are not included (GW-RPA) or they are included through 
the introduction of a XC hole around screening electrons 
only (GW-TDLDA) or around both the image-state elec- 
tron and screening electrons (GWF-TDLDA). We have 



reached the important conclusion that GW-RPA calcu- 
lations produce decay rates that are very close to GWF- 
TDLDA calculations, which are obtained with full inclu- 
sion of XC effects. With the use of either the free-electron 
mass or more realistic effective masses for all final states, 
decay rates of the n = 1 image state on Cu(lOO) are 
found to be smaller than those of the n = 1 image state 
on Cu(lll). The large bulk-state overlap on Cu(lll) is 
found to be approximately counterbalanced by the band 
gap extending on the Cu(lll) surface below the Fermi 
level, and differences between decay rates on the (100) 
and (111) surfaces of Cu are found to be mainly due to 
the large contribution, in the case of Cu(lll), from de- 
cay into the crystal-induced n — surface state located 
within the projected band gap. The results we have ob- 
tained for the n = 1 image-state lifetime on both sur- 
faces of Cu with use of realistic effective masses, which 
are ~ 20% below those obtained with mf = 1, are in ex- 
cellent agreement with experimentally determined decay 
times. 
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' The decay rate of Eq. 

" = 1 a.u. = 658 meV fs. 



by 



TABLE L GW-RPA decay rates, in linewidth units (meV), 
of the n = 1 ima ge s tate on Cu(lOO), as obtained from Eqs. 
(§) 5 (i) 1 (@) 1 s-iid ( p^ . Effective masses have been set equal to 
the free-electron mass (m/ = 1). All wave functions entering 
Eqs. (^) and ( |ll| ) have been chosen to be MP wave functions, 
and the final-state wave functions entering Eq. ^ have been 
taken to be either J, JG, or MP wave functions (see the text). 





Bulk 


Vacuum 


Inter. 


Total 


J 


21 


58 


-12 


67 


JG 


19.5 


13 


-11 


21.5 


MP 


24 


14 


-16 


22 



TABLE II. Decay rates, in linewidth units (meV), of the 
n = 1 image state on Cu(lOO), together with the experimen- 
tally determined decay rate of Refs. [l^ and p^ . All the wave 
functions entering Eqs. (^), (^ and ( |11[ ) have been chosen to 
be MP wave functions. Effective masses have been set equal 
to either the free-electron mass {rrif — 1) or to realistic values 
for all available final states (m/ 7^ 1). Five different models 
for the description of the screened intasaction have been used: 
the specular refiexioH. model (SRM)J2jthe model of Persson 
and Zaremba (PZ)J2j, and three self-consistent many-body 
approaches, GW-RPA, GW-TDLDA, and GWF-TDLDA. 



nif 


W 


Bulk 


Vacuum 


Inter. 


Total 


= 1 


SRM 


18 


3 


-4 


17 


= 1 


PZ 




55 




55 


= 1 


GW-RPA 


24 


14 


-16 


22 


= 1 


GW-TDLDA 


26.5 


14 


-16 


24.5 


= 1 


GWF-TDLDA 


23.5 


14 


-16 


21.5 


/I 


GW-RPA 


7 


11.5 


-1 


17.5 


/I 


GWF-TDLDA 


6.5 


11.5 


-1 


17 



Exp. 



16.5 



7 



TABLE III. As in Table II, but for Cu (111) and together 
with the deca^ rate experimentally determined for this sur- 
face of coppeiOat low temperature, T = 25K. Contributions 
to the linewidth from decay into bulk states lying below the 
bottom of the projected band gap, thereby excluding the con- 
tribution from decay into the n — intrinsic surface state, are 
displayed in parentheses. 



ruf 


W 


Bulk 


Vacuum 


Inter. 


Total 


Exp. 


= 1 


SRM 


46(34) 


12(1) 


-22(-5) 


36(30) 




= 1 


PZ 




57(2) 




57(2) 




= 1 


GW-RPA 


44(28) 


47(5) 


-54(-12) 


37(21) 




= 1 


GW-TDLDA 


43 


42 


-45 


40 




= 1 


GWr-TDLDA 


43.5 


47 


-54 


36.5 




/I 


GW-RPA 


32(24) 


34(5) 


-37(-12) 


29(17) 




/I 


GWr-TDLDA 


30.5 


35 


-38 


28.5 


30 



FIG. 1. Probability density of the n = 1 image state on 
Cu(lOO) (dotted line) and Cu(lll) (solid line), as obtained 
with use of the one-dimensional model potential of Ref. 
The crystal-edge {z = 0) has been chosen to be located half 
an interlayer spacing beyond the last atomic layer. Full circles 
represent the atomic positions. In the case of Cu(lll), the 
matching at the surface occurs at minimum amplitude. 

FIG. 2. Imaginary part of the n — 1 image-electron 
GW-RPA self-energy Im [— E(z,2';k|| = 0, i?o)] , versus z, in 
the vicinity of the (100) and (111) surfaces of Cu, as obtained 
from Eq. with use of MP wave functions both in Eq. 1^ 
and Eq. dill) and with all effective masses set equal to the 
free-electron mass, z' is fixed at —5 (top panel), (middle 
panel) and 5 A (bottom panel). 



TABLE IV. Calculated decay rates, in linewidth units 
(meV), and lifetimes, in femtoseconds (fs), of the n < 3 im- 
age states on Cu(lOO), together with the experimentally de- 
termined lifetimes of Refs. |l|and|l9[ In the case of n = 2 
and n = 3 contributions from decay into the lower lying im- 
age states are also displayed. The screened interaction has 
been evaluated within the GW-RPA. Model potential (MP) 
wave functions have been used for the description of all sin- 
gle-particle states entering Eqs. (^), (^, and (pi]). The effec- 
tive mass has been set equal to either the free-electron mass 
or to the realistic values described in the text. The result 
of introducing into Eq. (^) a realistic effective mass for all 
available final states is displayed in parentheses. 

Linewidths Lifetimes Experiment 



n = 1 

n = 2 

n = 2 to n 

n = 3 

n = 3 to n 

n = 3 to n : 



22(17.5) 
5(3.9) 
0.5(0.4) 
1.8(1.4) 
0.05(0.05) 
0.17(0.16) 



30(38) 
132(168) 



40 ±6 
120 ± 15 



367(480) 300 ± 20 



FIG. 3. Maximum of the imaginary part of the n = 1 im- 
age-electron GW-RPA self-energy Im [— E(2:, z'; ky = 0, -Eo)] , 
versus z' , in the vicinity of the (100) (dotted lines) and (111) 
(solid lines) surfaces of Cu, obtained as in Fig. 2. Vertical 
dotted and solid lines represent the atomic positions along the 
(100) and (111) directions, respectively. 



FIG. 4. As in Fig. 4, for the separate contributions to the 
magnitude of the maxiumum of Im [— E] of the n — 1 image 
electron on Cu(lll), according to wether transitions to bulk 
states (dotted line) or to the n = surface state (dashed line) 



FIG. 5. (a) Bulk, (b) vacuum, (c) interference, and (d) to- 
tal contributions to the GW-RPA damping rate coming 
from the decay into the various / available states lying be- 
low the projected band gap in Cu(lOO) (curves with circles) 
and Cu(lll) (curves with squares), as obtained with the fi- 
nal-state effective mass set equal to the free-electron mass 
(solid lines) and with use of realistic values of m/ (dashed 
lines). Model potential (MP) wave functions have been used 
for the description of all single-particle states entering Eqs. 

(i)' ^'^d (|ll|). Vertical dotted lines represent the Fermi 
level. 
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